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Abstract 

A family of potential-density pairs that represent spherical shells with 
finite thickness is obtained from the superposition of spheres with finite 
radii. Other families of shells with infinite thickness with a central hole 
are obtained by inversion transformations of spheres and of the finite 
shells. We also present a family of double shells with finite thickness. 
All potential-density pairs are analytical and can be stated in terms of 
elementary functions. For the above-mentioned structures, we study the 
circular orbits of test particles and their stability with respect to radial 
perturbations. All examples presented are found to be stable. A particular 
isotropic form of a metric in spherical coordinates is used to construct 
a General Relativistic version of the Newtonian families of spheres and 
shells. The matter of these structures is anisotropic, and the degree of 
anisotropy is a function of the radius. 

Key words: gravitation - galaxies; kinematics and dynamics. 

1 Introduction 

Analytical potential-density pairs with spherical symmetry have played an im- 
portant role in galactic dynamics. Simplified models may be used not only 
to test and to perform sophisticated numerical simulations, but also to help 
us gain insight into more complicated phenomena. Several spherical analytical 
potential-density pairs have been proposed as models for elliptical galaxies and 
bulges of disc galaxies [D [H |3] , for dark matter haloes (4j , as well as generaliza- 
tions of these models (SUHIII]. Shells constitute another spherical self-gravitating 
system that is important in issues like cosmology, gravitational collapse and su- 
pernovae (see, for instance, [8 and references therein). A common simplification 
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is to consider shells with infinitesimal thickness. Rein [5] constructed analytic 
solutions of the Vlasov-Poisson and Vlasov-Einstein systems representing, re- 
spectively, Newtonian and General Relativistic shells with finite thickness and 
a vacuum at the centre. A more detailed analysis of relativistic shells was done 
by [T^. Solutions for relativistic double and also multishells were found by [TT] 
using numerical methods. Other numerical solutions of finite Vlasov-Einstein 
shells were obtained by [12] . Recently, [13] proposed a simple Newtonian family 
of potential-density pairs of spherical shells with varying thicknesses and also a 
General Relativistic version of these models. 

In this work, we obtain families of potential-density pairs that represent 
spherical shells with finite thickness and shells with infinite thickness that have 
a central hole. We also build a family of double shells with finite thickness. For 
these structures, we study the rotation curves (circular orbits) of test particles 
and make a first stability analysis by considering the stability of circular orbits 
under small radial perturbations. We also consider a General Relativistic version 
of some of these Newtonian potential-density pairs in the same manner as in our 
previous work [T^. The paper is divided as follows. In Section [51 we present a 
family of spheres with finite radii that will be used as a basis for the construction 
of the spherical shells. In Section [31 the inversion theorem is applied on the 
spheres to obtain a first family of infinite shells. Then, in Section [H shells 
with finite thickness are constructed by superposing different members of the 
family of spheres. The inversion theorem is also used on these shells to obtain 
another family of infinite shells with a central hole. In Section[5l we show how to 
construct a family of double shells with finite thickness by superposing different 
members of finite shells, and we discuss a particular example. In Section [6] we 
use a particular isotropic form of a metric in spherical coordinates to construct 
a simple General Relativistic version of the Newtonian spheres and shells. We 
study in some detail two examples of spherical structures. Finally, in Section [7] 
we summarize our results. 

2 A family of finite spheres 

We begin by considering a family of spheres with finite radius and following 
mass-density distribution, 



where m = 1,2,.... This form of mass distribution has been chosen because 
it is a monotononically decreasing function of the radius and the corresponding 
potential can be expressed in terms of elementary functions, as will be seen 
below. Furthermore, models of thin finite discs with surface densities in the 
form of ll]), where r now represents the cylindrical radius, are well known. 
The disc member with m = 1 is the uniformly rotating disc of Kalnajs [14] . 
Gonzalez & Reina [H] constructed a family of generalized Kalnajs discs with 
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surface density similar to (HJ; however, it is interesting to note that a family of 
discs with a similar density distribution has been used much earlier by Morgan 
& Morgan [TB] in their study of General Relativistic static thin discs. 
The total mass Mm of the sphere is 



Mr, 



47r / r pjndr 



PcTT^a^{2m~ 1)!! 
2'"(m + 1)! 



(2) 



For convenience, we consider spheres with the same mass M — Mm- Then, 
equation ([T]) can be rewritten as 



M2"(to + 1)! 
7r2a3(2m- 1)!! 



2 ^ m-1/2 



(3) 



The potential of the sphere with m = 1 follows from direct integration of the 
Poisson equation in spherical coordinates: 



j__d 
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,d$\ 
dr J 



= 4:TtGp, 



and is given by 
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The potentials of the spheres with m > 1 can be found by using the recurrence 
relation 

'i>™+i - I «^'"+^<f™da. (6) 



For reference, the potentials of the members with m — 2,3,4: read as 
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arcsin (^-) - Va^ - (48r^ - 200r'*a2 + 326r2a'' - 279a'^) 

(8) 



315ai 



■ arcsin 



(V^ + Va2~72 (I28r* - GSGr^a^ + ISeSr-^a"* 

-WQOr^a^ + 965a^)] . (9) 

A common dynamical quantity used to characterize spherical distributions 
of matter is the circular velocity Vc (or the rotation curve) of test particles. 
We will mosty be interested in the circular orbits of particles inside the matter 
distributions. We also investigate the stability of these orbits under small radial 
perturbations. This is characterized by the epicyclic frequency k, where stable 
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Figure 1: (a) The mass density pm — Pm/{M/a^), equation ([3]), as func- 
tion of r/a, for four members of spheres, (b) The circular velocity Wc(m) = 
Vc{m)/{GM/ay^^, equations (j46l) - (|49|) . and (c) the epicyclic frequency Rm = 
Km/{GM/a^y^'^, equations (|50l) " (l53|) . as functions of r/a, for the same mem- 
bers of spheres. 



circular orbits correspond to > 0. For spherical potentials, the circular 
velocity and epicyclic frequency are calculated through the relations [17] 

Vc = ^/r¥~^, and n'^ = ^ + ■ (10) 

Expressions for the circular velocity Wc(in) ^-nd epicyclic frequency Km for the 
spheres with m = 1, 2, 3, 4 are listed in Appendix 1X1 Outside the spheres {r > a), 
we have = [GM/rf/^ and k = [GM/r^Y'^. 

In Fig. [TJa), we plot some curves of the dimensionless mass density pm = 
Pm/{M/a^), equation ([3]), as function of r/a, for four members of the family of 
spheres. In Figs [Ijb) and (c) we show, respectively, some curves of the dimen- 
sionless circular velocity Vc{m) = Vc(m)/iGM/ay/^, equations (|46l) -(|49 | . and 
curves of the dimensionless epicyclic frequency Km — K„i/(GM/a^)^/^, equa- 
tions ((50)) - 153]) . as functions of r/a, for four members of the family of spheres. 
Unlike the disc case, we do not have a linear rotation profile for the sphere with 
m = 1. We also note that the circular orbits for the four first members are 
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stable. 



3 Inverted spheres 

A new family of spherical potential density-pairs can be obtained by making an 
inversion (or a Kelvin transformation) [T51[TH] on the spheres with density ([3]). 
In spherical coordinates, the inversion theorem states that if a potential-density 
pair p(r), $(r) is a solution to the Poisson equation, then the pair 

P{^nv:) = (^) P , (11) 

is also a solution to the Poisson equation. The inversion theorem applied to the 
density ([3]) results in 

M2™a2(m + 1)!/ a^x^-i/s 

This mass distribution can be interpreted as a family of spherical shells with 
a central hole of radius a. The shells extend to infinity, although the density 
decays quite fast, Pm{inv.) oc The maximum of density occurs at r/a = 

yj {2m + 4)/5. The mass Mm of each shell is 

M,n = 47r / r^p„(,™.)dr = ^^..^ -ttj. 14) 

^ ' 7r(2m + l)(2m — 1)!! 

For convenience, we consider shells with the same mass M = Mm a-nd rewrite 
([T3D as 

(2m + l)Xa2/ 



Pm(^nv.) = \'^) ' - (1^) 

The expressions for the potentials of the shells with m = 1, 2, 3, 4 are given 
in Appendix [BJ From them result simple equations for the rotation curves 
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Vc{m){2nv.) and epicyclic frequencies namely 
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which suggest the general relations 
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Some curves of the dimensionless mass density Pm{inv.) = Pm{inv .) / {-^ / '^^) j 
equation (fTS]) . as function of r/a, are shown in Fig. [^Ja). The shell with m = 1 
exibits a quite sharp density profile. In Figs[2Ib) and (c), we plot, respectively, 
curves of the dimensionless circular velocity Vc(m)(inv.) = Vc{m){inv .) / {GM / a)^ , 
equations (|16p and P?|) . and the dimensionless epicyclic frequency ) = 

i^m{inv.) / (GM / a^y^^ , equations (fT8 l) - (f2T1) . as functions of r/a. The four mem- 
bers of the shells shown are stable with respect to radial perturbations of circular 
orbits. From we see that all members are stable. According to Newton's 
shell theorem, the potential inside the hole (0 < r < a) is constant and thus the 
circular velocity must be zero. 



4 Families of spherical shells 

Now we use the spheres discussed in Section [2] to construct potential-density 
pairs that represent shells of matter with finite thickness. The method is similar 
to the used by [20] to obtain rings by superposing the Morgan & Morgan discs 
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Figure 2: (a) The mass density Pm(mu.) = Pm{inv.)l{M^/o!^)j equation 
as function of r/a, for four members of inverted spheres, (b) The circular 
velocity Vc{ra){inv.) = Vc{m){inv.) I {GM / o)^/'^ , cquatious (HH) and and (c) 

the epicyclic frequency K,m{inv.) — '^m{inv.)/{GM/a'^y^^, equations P^ - (PT|) . 
as functions of r/a, for the same members of inverted spheres. 
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[IB] . We take the density distribution ([U and make the fohowing sum: 



(m,n) 



k=0 



2X^+1/2 n / „2\"-fe / ^2\™+l/2 2« 



(24) 



where = n!/[fc!(n — fc)!], m — 0,1,... and n — 1,2, This new den- 
sity distribution vanishes on r = and r = a, so it represents a spheri- 
cal shell with thickness a. Its mass density has a maximum value at r/a = 
y/{2n)/{2m + 2n+l). The total mass M^™^") of each sheh is 

M(-") . 4. /" r^pi--Ur = P^^^^^j^m + l)!!(2n + 1)!! _ 

Jo 2™+"+i(TO + n + 2)! ^ 

The potentials associated with the mass density (1^^ are found by a similar 
sum, i.e. the potential of the shell with m = 0,n = 1 is = — $3. The 

explicit expressions for the potentials of the shells with rn ~ 0,1 and n — 1,2 
are given in Appendix ICl 

In Fig. [Sja), we display curves of the dimensionless mass density p^™'"' = 
equation (|24|) . as function of r/a, for four members of spherical shells. 
For small values of r/a, the density of the shell grows as (r/a)2("+i); therefore, 
shells with larger values of n also have larger holes. Curves of the dimensionless 
circular velocity Wc™'"^ = Wc™'"V(GPca^)^^^, equations (|62l) - ([65|) . and of the 
dimensionless epicyclic frequency /tt™'") = K^™'"V(GPc)"^^^, equations (|66| - 
are displayed in Figs[3Jb) and (c), respectively. For small values of r/a, 
the rotation curves are proportional to (r/a)"+^ and thus also grow more slowly 
as n increases. We also see that the members of shells shown in the figure are 
stable. 

It is also possible to invert the shells. A Kelvin transformation applied to 
the density ([^ results in 



a^\ ' /a\2"+5 



P&:!-Pc(l-^j (^j ; r>«, (26) 

which represents a shell with a central hole of radius a. If n = and m — ^ m — 1 
we recover the family of inverted shells discussed in Section [3l They also extend 
to infinity with an even larger decay of density, p^™!,^) oc 1/ r^^^^ . The maximum 

of density occurs at r/a — \J (2m + 2n -t- 6)/(2n -|- 5). The total mass of each 
shell is 

^("'^■) (2m + 2n-K3)!! ^ ' 

Curves of the dimensionless mass density equation (j26p . as function 

of r/a, are displayed in Fig. lU^a) for four members of inverted shells. Figs 



8 




Figure 3: (a) The mass density p^™'"^) = p'^™-''^'^ / p^^ equation (|24l) . as function 
of r / a, for four members of spherical shells, (b) The circular velocity = 
v^'^'' / {GpcO^Y^'^ , equations (|62)) ~ (|65|) . and (c) the epicyclic frequency k^™-") = 
((7^^)1/2^ equations ([66|) - (|69]) . as functions of r/a, for the same members 
of spherical shells. 
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Figure 4: (a) The mass density = equation ((261), as func- 

tion of r/a, for four members of inverted shells, (b) The circular velocity 
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[DJb) and (c) show, respectively, curves of the dimensionless circular velocity 
_(m,n) gq^a,tions (|74p -(|77 p . and of the dimensionless epicyclic frequency k^'"'"-' 



equations ([75|) - ([5T|) . for the same members of inverted shells. Once again, these 
members are stable with respect to radial perturbations of circular orbits. 

5 A family of double shells 

The members of the family of shells discussed in Section 0] may be combined 
to generate double shells with finite thickness. To achieve this, we superpose 
members of (j24p in the following way: 

(28) 

where 6 > 1. This density distribution may be interpreted as two concentric 
spherical shells with a gap located at r/a = 1/b. We take as an example the 
member with m — 0,n — 1. The potential as well as the expressions of the 
rotation curve and epicyclic frequency for this example are given in Appendix 

m 

Fig. El^a) shows curves of the dimensionless mass density P(2)^'' — P'(^'if' I Pc, 
equation (|28l) . as function oi r/a for some values of b. Curves of the circular 
velocity v'f^2) ~ '^c(2)V(GPca^)^/^, equation (l83l) . and of the epicyclic frequency 
_(m,n) _ ^(™'")^^(];p^^i/2^ equation (f84|) . are plotted in Figs^b) and (c), re- 
spectively. Unlike the structures discussed so far, the rotation curves for this 
example of double shell show a point of minimum. They also show no sign of 
instability. 



6 General Relativistic Spheres and Shells 

Now we consider an extension to General Relativity of the Newtonian spherical 
potential-density pairs discussed so far. We choose a particular metric in an 
isotropic form in spherical coordinates (t, r, 9, ip), which was also used in a recent 
work on another relativistic model of shells [13], 

ds2 = ' c^df -{1 + ft (dr^ + r^de^ + sin^ Odcp^) , (29) 

where / = /(r). The Schwarzschild solution is given by metric ([^^ if / = 
GAf/(2c^r), and the relation between the function /(r) and the Newtonian 
potential $(r) is 

f-~l. (30) 
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Figure 5: (a) The mass density p^2) = P(2) I Pc^ equation as function 

of r/a, for a double shell, (b) The circular velocity v^^(^2) ~ ^c(2)V(GPca^)^^^, 
equation (|83l) . and (c) the epicyclic frequency — k|^'"V(G'/9c)^^^i equa- 

tion (|84p . as functions of r/a, for the same double shell. 
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By using the Einstein equations, we find the following expressions for the non- 
zero components of the energy-momentum tensor T^^ |13) 



* 27rG(l + ff dr V dr 



d/ // , d/\ 



rpip 



2i:G (1 -f ff (1 - /) dr V J- dr / 



,d2/ , /d/ 



4^G(1 + /)^(1-/) 



dr 



(31) 
(32) 

(33) 



The energy density reads as e = Tl/c^ and the pressures or tensions along a 
direction k are given by Pk = —T^- The 'effective Newtonian density' pn = 



e + Pr/c + Pe/c + Pip/c can be cast as 



1 d 



Pn = 



27rG{l + fy (1 - /) dr 



>d/ 
dr 



(34) 



We will study in some detail two simple examples of General Relativistic 
spherical structures with finite extent: the sphere with potential ([5]) and the shell 
with potential ([58l) . Using equations ([5|) and (|30l) - (|34l) . we find the following 
expressions for the components of the energy-momentum tensor of the sphere, 



^2a3(i + /i)= 
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(38) 



(39) 



and the dimensionless variables and parameters are f — r/vs, a — a/vs, ei = 
ei/(M/r3), PNl = pjvi/(M/r3), P^, = Pki/{,M^/rl) and r, = GAf/(2c2). 

In order to represent physically meaningful matter distributions, the com- 
ponents of the energy- momentum tensor should satisfy the energy conditions. 
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Figure 6: (a) The energy density ei = ei/{M/r^g) (solid curve), equation 
(|35|) . and the effective Newtonian density pNi = pNi/{M/r^) (dashed curve), 
equation (pG)) . for the sphere with radius a = 7. (b) The radial pressure 
Pri = Pri/iMc^ /r^) (soUd curve), equation (|37l) . and the azimuthal/polar 
pressure Pgi = Pbi/{Mc^ /r^) (dashed curve), equation ([55)1 . for the sphere 
with radius a = 7. 



The strong energy condition states that > 0, whereas the weak energy 
condition imposes the condition e > 0. The dominant energy condition requires 
|-Pfc/e| < c^- Equations ([55]) and show that the strong and weak energy con- 
ditions are satisfied if a > 16/(37r) « 1.70. In Fig.[6l^a), we plot the curves of the 
energy density ei (solid curve) and the effective Newtonian density pNi (dashed 
curve), as functions of f = r/r^, for the sphere with radius a — 7, and in Fig. 
[6fb), we display the radial pressure Pri (solid curve) and the azimuthal/polar 
pressure Pgi, as functions of f = r/r^, for the same sphere. Unfortunately, we 
found that the energy density is not always a monotone decreasing function of 
the radius. This happens if a ^ 6.8. For a ^ 6.8, the dominant energy condition 
is also satisfied. The radial and azimuthal pressures also decrease monotonically 
and both vanish on the sphere's surface. In this example, we have an anisotropic 
sphere with larger azimuthal/polar than radial pressures. 

In the case of the shell with potential (|58p . we obtain the following expres- 
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sions for the components of the energy-momentum tensor: 
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and the dimensionless variables and parameters are now f ~ t/tc 
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From equations (|40|) and (|4T|) . we have that the strong and weak energy 
conditions are satisfied when a < \/T5/2 ss 1.94. The dominant energy condition 
holds if 



a < 



5V67r + 15 



1.29. 



(45) 



2(27r + 5) 

Fig-EIa) shows curves of the energy density e^O'^) (solid curve) and the effective 
Newtonian density Pn'^^ (dashed curve), as functions of r/rc, for the shell with 
parameter a = 1.2. In Fig. El^b), we show curves of the radial pressure Pr^'^^ 
(solid curve) and the azimuthal/polar pressure P^°'^^ (dashed curve), as func- 
tions of r/vc, for the same shell. The pressures, as well as the densities, vanish 
at the centre and and at the outer surface of the shell. This example shows 
a highly anisotropic relativistic shell with larger azimuthal/polar than radial 
pressures. We note that the analytical and numerical solutions for relativistic 
finite shells found by [TU] and [H] also have anisotropic pressures. 



7 Discussion 

In this work, we constructed families of potential-density pairs that represent 
spherical shells with finite thickness and shells with infinite thickness, but with 
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Figure 7: (a) The energy density e^"'-'^^ = e^^'^^/pc (solid curve), equation 
pO|) . and the effective Newtonian density p^'^'' = Pn'^'^/Pc (dashed curve), 
equation (|4T|) . for the shell with a — 1.2. (b) The radial pressure pj^'^-* = 
Pr^''^^ /(pcC^) (solid curve), equation P^ . and the azimuthal/polar pressure 
Pjj'^'^^ — Pg^'^^ / [p^c^) (dashed curve), equation (|43| . for the shell with a — 1.2. 



a hole at the centre. Shells with finite thickness were obtained by the superpo- 
sition of spheres with finite radii, and shells with infinite extend by using the 
inversion theorem on spheres and finite shells. Furthermore, shells with finite 
thickness may be superposed to give a family of double shells with finite thick- 
ness. All the potential-density pairs discussed in this work can be expressed in 
terms of elementary functions. We also studied the rotation curves for these 
self-gravitating structures and made a stability analysis based on radial pertur- 
bations of circular orbits of individual particles. In all the examples discussed, 
we found stable structures. 

A General Relativistic version of the Newtonian spherical potential-density 
pairs was obtained by using a particular isotropic form of a metric in spherical 
coordinates. We studied in some detail an example of a relativistic sphere and of 
a shell. In both cases, the matter presents equal azimuthal and polar pressures, 
but different from the radial pressure, and this anisotropy is a function of radius. 
It is possible to choose parameters such that all energy conditions are satisfied. 
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A Rotation curves and epicyclic frequency for 
spheres 

The expressions for the circular velocity i'c(m) S'Hd epicyclic frequency Km for 
the spheres with m = 1, 2, 3, 4 read as 



"c(l) 



Vc{2) 
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"c(3) 



2GM 
157ra^ 



2GM 
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15a8 



— arcsin ( - ) + \J o? — r"^ (2r^ — 
r Va/ ^ 



1/2 



— arcsin (-) - ^ o? - (8r^ - Mr^a^ + 3a^) 



1/2 



(46) 
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arcsin (-) - - r"^ (384r^ - 1488r^a2 



+2104r4a4 - 1210r2a'^ + lOSa^)] , 



Ki 



2GM 



— arcsin 
r Va 



in + ^/a? - r2 (lOr^ - a^) 



1/2 



K2 = 



«3 



2GM 
37rr2a6 

2GAf 
157rr^a® 



— arcsin ( -) - ^a^ - r"^ (SSr^ - 62r'a2 + 3a*) 



1/2 



(49) 
(50) 

, (51) 



15a** 



arcsin + s/ - (432r^ - 904r*a2 + 502r2a'' - 150^*) 

(52) 



/t4 



2GM 
lOSTrr^aio 



105a 



10 



arcsin (0 - Va^ - (4224r^ - 13008r^a^ 



+13624r*a'' - 5050r^a^ + 105a^)] 



8M1/2 



(53) 



18 



B Potentials of inverted spheres 

The expressions for the potentials ^rn{inv.) of the sheUs (|15p with m — 1,2, 3, 4, 
obtained by inversion of spheres, are 
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C Potentials, rotation curves and epicyclic fre- 
quencies for spherical shells 

The expressions for the potentials of the shells (p4)) with m — 0,1 and 

n — 1,2 are 
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Wc also give the expressions for the circular velocity u^'"'"-' and epicyclic fre- 
quency K^™'") for these members: 
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D Potentials, rotation curves and epicyclic fre- 
quencies for inverted spherical shells 

The expressions for the potentials of the inverted shehs (p6)) with m = 0, 1 



and n = 1,2 are 
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The circular velocity v^J^^J^ and epicyclic frequency for these inverted 

shells read as 
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E Potential, rotation curve and epicyclic fre- 
quency for a double shell 

The potential ^"(2)^^ double shell discussed in Section [5] reads as 
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The corresponding expressions of the circular velocity v^^(^2) ^"^^ epicyclic fre- 

(0,1) 

quency ^^2) ^-re 



^(0,1) _ , / GpcTT j 15a 



10 



(76^ - 206^ + 16) arcsin (^-) 



^=(2) V 960a7 

- r2 [384r*^&'^ - 4:Sr^a^b^ {b^ + 20) - Sr'^a'^ (76'* - 205^ - 8O) 
-lOr^a^ (754 _ 2062 + 16) - 15a« (76^ - 2062 + 16)] j'^' , (83) 



(0.1) _ / GpcTT 1 15a 
"(2) " V 960r2a 



10 



^76^^ - 2062 + 16) arcsin (^-) 



+ Va2 - r2 [4224r^6^ - 48r'^a262 (62 + I8O) - 8r^a^ {ib^ - 20b^ - 560) 

-10r2a6 (7^4 _ 20^2 + ig) _ 15^8 (76^ - 2062 + 16)] j'^' . (84) 



22 



